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Abstract. A QCD analysis of the world data on inclusive polarized deep inelastic scattering of 
leptons on nucleons is presented in leading and next-to-leading order. New parameterizations are 
derived for the quark and gluon distributions and the value of a, {Mz) is determined. Emphasis is put 
CO . on the derivation of fully correlated error bands for these distributions which are directly applicable 

to determine experimental errors of other polarized observables. The impact of the variation of both 
the renormalization and factorization scales on the value of a^ is studied. Finally a factorization- 
?*■ . scheme invariant QCD analysis based on the observables g\{x,Q^) and dg\{x,Q^) / d\og{Q^) is 

performed in next-to-leading order, which is compared to the standard analysis. 
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5 INTRODUCTION 

"Th ! The remarkable growth of experimental data on inclusive polarized deep inelastic scat- 

P^l tering of leptons on nucleons over the last years [1-9] allows to perform refined QCD 

p ■ analyses of polarized structure functions in order to reveal the spin-dependent partonic 

(D ! structure of the nucleon. A number of such analyses has already been worked out. The 

'^ ' most recent ones are [10-13] \ In this talk results from a new QCD analysis in leading 

(LO) and next-to-leading (NLO) order [14] ^ are presented. New parameterizations of 
^ . the polarized quark and gluon distributions are derived including the parameterizations 

of fully correlated lo error bands for these distributions, which are directly applicable 
to calculate errors of other polarized observables. Furthermore the value of as{Mz) is 
determined. Finally and for the first time a factorization-scheme independent QCD evo- 
lution based on the observables gi{x,Q^) and dgi{x,Q^)/d\og(Q^) in next-to-leading 
order is performed. 



FORMALISM 

In LO the polarized structure function ^i (x, Q^) is expressed as the sum of the polarized 
quark distributions A^,(jc, Q^) weighted by the square of the quark charges. In NLO the 
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^ For a more complete list of references see the references therein and Ref. [14] 
All details of the analysis are given in Ref. [14]. 



expression for gi{x, Q^) involves the polarized singlet AE(x, Q^), the gluon AG(x, Q^ 
and the non-singlet /^q'^^{x, Q^) distributions and reads 



gx{x,Q^) = \ 



(1) 



— E"- 1 4] [SQ ® AE + 5Cg ® AG] + 5Cw5 ® A^^^ 

where n/ is the number of active quark flavors and ei is the quark charge. The symbol 
denotes the Mellin convolution w.r.t. x of the polarized parton densities /^qi{x, Q^) with 
the corresponding polarized Wilson coefficient functions 5C/(x, as(2^)). The polarized 
singlet and non-singlet distributions are certain combinations of the polarized quark 
distributions /^qi{x, Q^). 

The evolution equations used to evolve the parton densities to different (^ values 
contain the polarized splitting functions AP,j(x, av(2^)). Both the polarized Wilson 
coefficient [15] and the polarized splitting functions [16] are known in the MS scheme 
up to order O (a^) . 



METHOD 

The shape chosen for the parameterization of the polarized parton distributions at the 
input scale of (^ = 4.0 GeV^ is : 

xAqiix, QI) = ri,-A,-.)c"'(l -xf'{\+yiX + PiX^). (2) 

The normalization constant A, is chosen such that r|; is the first moment of A^,(jc, 2q). 
The densities to be fitted are Am^^, AJy, A^, and AG. 

Assuming SU{3) flavor symmetry the first moments of Amv and AJy are determined by 
the SU{3) parameters F and D measured in neutron and hyperon |3-decays and can be 
fixed to r[u^ = 0.926 and Tj^^, = —0.341. In addition we assume a flavor symmetric sea, 
i.e. only one general sea distribution Aq{x, Q^) is required. No assumptions are made 
concerning positivity and helicity retention. Given the present accuracy of the data we 
set a number of parameters to zero, namely p„^, = p^/^, = 0, y^ = p^ = 0, and Yg = Pg = 0. 
This choice reduces the number of parameters to be fitted for each parton distribution 
to three. In addition the parameter A^co was determined. The relative normalizations 
of the different data sets were fitted and then fixed. Doing so part of the experimental 
systematics was taken into account. 



RESULTS 

The results reported here are based on 433 data points of asymmetry data, i.e. gi/Fi 
or Ai, above Q^ = 1.0 GeV^, the world statistics published so far. The QCD fits are 
performed on ^i which is evaluated from the asymmetry data using parameterizations for 



Note that: A^ + A^ = Aq^, + 2Aq. 



the unpolarized structure functions F2 [17] and 7? [18]. We realized that the 4 parameters 
liiy^ Jdy^ t>q, and bo had to be fixed in addition at their values at %^-j^ since the data do 
not constrain these parameters well enough. Only fits with a positive definite covariance 
matrix were accepted in order to be able to calculate the fully correlated lo error bands. 
The NLO polarized parton densities at the input scale are presented in Fig. 1 . 
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FIGURE 1. 



Polarized parton distribution at the input scale (A 



4.0 GeY (solid line) compared to 



results obtained by GRSV (dashed-dotted line) [12] and AAC (dashed line) [10]. The shaded areas 
represent the fully correlated la error bands calculated by Gaussian error propagation, Ref. [14]. 

While the quality of the data is sufficient to determine Awy and AJv with good 
accuracy, AG and A^ have much broader error bands. This is essentially due to the lack 
of data at low x. The agreement with the results of the analyses of Refs. [10] and [12] is 
satisfactory within the error bands. The measured structure function g^ is well described 
both as function of x and of Q^ . The derived parton distributions and its error bands have 
been evolved to Q^ values up to 10,000 GeV^ . As an example the evolution of AG is 
shown in Fig. 2. One observes that even within the error band AG stays positive up to 
the highest Q^ value. It should be mentioned that t^q develops a trend to change sign and 
becomes slightly positive towards higher Q^ values and for xJc;0. 1 within the errors. 
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FIGURE 2. The polarized parton distribution AG evolved up to Q^ values up to Q^ ~ 10,000 Gev^ 
(solid line) compared to results obtained by GRSV (dashed-dotted line) [12] and AAC (dashed line) [10]. 
The shaded areas represent the fully correlated la error bands calculated by Gaussian error propagation, 
Ref. [14]. 



In determining qLs{M^ the parameter AgcD was fitted. The impact of the variation of 
both the renormalization and factorization scales on the value of a, was studied. The 
following value for Aqcd was obtained 



A, 



(4) 



+ 117 



^^6D = 241±58(fit) :;4(fac, _^^ 
which results into a value of 



fren) MeV, 



+0.004 ,^ , +0.005 



+0.008 



a,(M2) =0.114 _q;qq3 (fit) _q;qq4 (fac) _^_^^^ (ren 



This value of a.5(M|) is compatible within the errors with the world average of 0. 1 18 ± 
0.002 [19] and with values from other QCD analyses [20], although the central value 
tends to be lower, as also in Ref . [20b] . 

Finally a factorization-scheme invariant QCD analysis based on the observables 
gi{x,Q^) and dgi{x,Q^)/dlog{Q^) in next-to-leading order was performed. The cor- 
responding evolution equations have been worked out in Ref. [21] ^. Such an analysis 
has the advantage of direct control over the input since it comes from measured quan- 
tities. The only parameter to be determined is Aqcd- Unfortunately, the present data do 
not yet allow to determine the slope dgi (x, 2^)/5log(2^) as an input density from mea- 
surements of ^1, but it is derived here from the fit result for gi described above. The 
evolution of the so determined slope is shown in Fig. 3. 
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FIGURE 3. The evolution of dg\{x,Q^)/dt (singlet contribution) with t = -2/poln(a,(e^)/a,(go))- 
The slope of ^i was determined from a fitted gi, see text. 



The same case has already been considered in Ref. [22]. 



A downward shift of 12 MeV in Agco was found yielding a similar result for av(M|) as 
obtained in the standard analysis. 



CONCLUSIONS 

An LO and NLO QCD Analysis of the current World-Data on Polarized Structure Func- 
tions was performed. New parameterizations of the polarized parton densities including 
their errors were derived. They are available via a fast FORTRAN code for the range: 
I <Q^ < 10^ GeV^ and 10""^ <x<l. The value determined for a,(M|) is compatible 
with the world average, although the central value obtained is lower. First steps in a 
factorization-scheme invariant QCD evolution based on the structure function ^i {x, Q^) 
and dgi (x, 2^)/3log Q^ were performed yielding similar results for a^(M|). This latter 
analysis is a very promising way to proceed in the future, since it allows to extract Aqcd 
fixing all the input distributions by direct measurements. 
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